Accepted by the Astrophysical Journal 

Optical Spectroscopy of GRO J1655-40 

Roberto Soria^, Kinwah Wu^, Richard W. Hunstead^ 

ABSTRACT 

We have obtained optical spectra of the soft X-ray transient GRO J1655— 40 during 
different X-ray spectral states (quiescence, high-soft and hard outburst) between 1994 
August and 1997 June. Characteristic features observed during the 1996-97 high-soft 
state were: a) broad absorption lines at Ha and H/3, probably formed in the inner disk; 
b) double-peaked He II A4686 emission lines, formed in a temperature-inversion layer on 
the disk surface, created by the soft X-ray irradiation; c) double-peaked Ha emission, 
with a strength associated with the hard X-ray flux, suggesting that it was probably 
emitted from deeper layers than He II A4686. The He II A4686 line profile appeared 
approximately symmetric, as we would expect from a disk surface with an axisymmetric 
emissivity function. The Balmer emission, on the other hand, appeared to come only 
from a double-armed region on the disk, possibly the locations of tidal density waves or 
spiral shocks. The observed rotational velocities of all the double-peaked lines suggest 
that the disk was extended slightly beyond its tidal radius. 

Three classes of lines were identified in the spectra taken in 1994 August-September, 
during a period of low X-ray activity between two strong X-ray fiares: broad absorp- 
tion, broad (fiat-topped) emission and narrow emission. We have found that the narrow 
(single-peaked or double-peaked) emission lines cannot be explained by a conventional 
thin accretion disk model. We propose that the system was in a transient state, in 
which the accretion disk might have had an extended optically thin cocoon and signif- 
icant matter outflow, which would also explain the systematic blue-shift of the narrow 
emission lines and the flat-topped profiles of the broad emission lines. 

After the onset of a hard X-ray fiare the disk signatures disappeared, and strong 
single-peaked Ha and Paschen emission was detected, suggesting that the cocoon be- 
came opaque to optical radiation. High-ionization lines disappeared or weakened. Two 
weeks after the end of the flare, the cocoon appeared to be once again optically thin. 
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1. INTRODUCTION 

Among the low-mass X-ray binaries, soft X-ray transients (SXTs) offer a favorable opportunity 
of detecting stellar-mass black-hole candidates (BHCs). Measurements of the binary periods and 
the orbital velocities of the companion stars in seven SXTs yield mass functions ^ 3Mq for their 
compact stars (McClintock 1998, Filippenko et al. 1999). As these masses are higher than the 
upper limit to the mass of neutron stars for known equations of state of nucleon matter, the 
compact stars in SXTs are all BHCs. In other two systems, the mass functions are < 3Mq but the 
inclination angles are thought to be sufficiently low that the compact objects can also be classified 
as BHCs (McClintock 1998). 

SXTs show transitions between different spectral states. In a simple classification, we can 
identify an "off" (or "quiescent" ) state when the source is inactive, and a "high-soft state" and a 
"low-hard state" defined by their X-ray spectral properties, when the source is X-ray active. Often 
an increase of the optical brightness and radio flares accompany the X-ray outbursts. (See e.g. 
Tanaka & Lewin 1995 for a review of the properties of SXTs.) 

GRO J1655-40 is an SXT discovered by BATSE on 1994 July 27 (Zhang et al. 1994). Between 
1994 and 1995 it underwent several X-ray outbursts (see Tavani et al. 1996), after which it retreated 
to a quiescent state. In late April 1996, another outburst started, and it lasted more than one year. 
The 2-12 keV X-ray light curve of GRO J1655— 40 during the 1996 outburst is shown in Figure 1 
(from the quick-look results provided by the RXTE/ASM team); in the same figure we also show 
the hard (20 - 100 keV) X-ray flux detected by BATSE over the same period of time. 

The binary period of GRO J1655-40 is P = 2.62157 ± 0.00015 d (Orosz & Bailyn 1997). The 
time of inferior conjunction of the secondary star is HJD 2449838.4209 ± 0.0055 (Orosz Sz Bailyn 
1997). An alternative ephemeris, given by van der Hooft et al. (1998), is HJD 2449838.4198(52) -|- 
2.62168(14) X N. 

The mass function determined from the kinematics of the system by Orosz Sz Bailyn (1997) 
implies a mass Mi = (7.0 ± 0.2) Mq for the compact object. By taking into account the effect 
of X-ray irradiation of the secondary star, Phillips, Shahbaz & Podsiadlowski (1999) derived lower 
values for the mass of the primary: they obtained 4.1 < Mi < 6.6Mq (90 per cent confidence 
level). By using only data taken during an X-ray quiescent state, Shahbaz et al. (1999) obtained 
5.5 < Ml < 7.9Mq (95 per cent confidence level). In any case, the mass of the compact object 
makes GRO J1655— 40 a candidate black hole binary. The major orbital parameters of the system 
are listed in Table 1. 

During the quiescent states the system is optically faint, with V ^ 17 mag (Bailyn et al. 1995b), 
but it can brighten substantially when an X-ray outburst occurs. In August 1994, it reached 
V = 14 A (Bailyn et al. 1995a), and in May 1996, V = 15.4 (Home et al. 1996). 

Since its discovery, we have observed GRO J1655— 40 spectroscopically on four occasions: 1994 
August - September, 1996 April, 1996 June and 1997 June. The observations covered periods in 
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which the system was in a pre-outburst quiescence, at the onset of a hard X-ray flare and in a 
high-soft state. In a previous paper (Soria et al. 1998) we presented the orbital-phase-dependent 
velocity shifts of the He II A4686 and N III AA4641,4642 emission lines. 

In this paper we extend our study to include the Balmer lines and the orbital variations in their 
line profiles. Technical details of our observations are outlined in §2. In §3 we present the results of 
the spectroscopic observations we conducted in 1996 April, and in §4 we present those we obtained 
in 1996 June and 1997 June. In §5 we discuss the profiles of the Balmer lines observed in the 1996 
- 1997 soft X-ray outburst, and we attempt to locate the emission regions. We also discuss some 
physical mechanisms responsible for line emission and absorption in the disk. In §6 we present the 
results of the observations we carried out during the 1994 hard X-ray outburst. We outline the 
general spectral features, we discuss the origin, intensity and profile of the main absorption and 
emission lines, and we compare the two sets of spectra from 1994 and 1996. In particular, we show 
the evolution of the optical spectrum before, during and after a major hard X-ray flare. 



2. OBSERVATIONS 

GRO J1655— 40 was observed by one of us (RWH) as a target of opportunity on each night 
from 1994 August 30 to 1994 September 4, with the RGO spectrograph and Tektronix Ik x Ik 
thinned CCD on the 3.9 m Anglo- Australian Telescope (AAT). Spectra were obtained in two 
regions, 6278 — 6825 A, centered on the Ha line, and 4432 — 5051 A, covering the N III, He II and 
H/? lines. Gratings with 1200 grooves/mm were used, with the blaze direction oriented towards the 
25 cm camera, giving a resolution of 1.3 A FWHM. On September 6, we obtained simultaneously 
a blue (3925 — 5500 A) spectrum with the RGO spectrograph and 600 grooves/mm gratings, at a 
resolution of 2.5 A FWHM, and a red (5500 - 11000 A) spectrum with the FORS spectrograph, 
via a dichroic beam splitter, at a resolution of 20 A FWHM. The latter observations were kindly 
obtained for us by Paul Francis. 

In 1996 April 20 - 21 we carried out a scheduled observation with the Double Beam Spectro- 
graph (DBS) on the ANU 2.3 m Telescope at Siding Spring Observatory. The detectors on the 
two arms of the spectrograph were SITe 1752x532 CCDs. Gratings with 300 grooves/mm were 
used for the blue (3600 — 5700 A) and the red (5700 — 9300 A) bands, and low-resolution spectra 
(resolution = 4.8 A FWHM) were obtained. More extensive observations were carried out in 1996 
June 8-12 and June 17, and 1997 June 14 - 15, with the DBS on the ANU 2.3 m Telescope. 1200 
grooves/mm gratings were used for both the blue (4150 — 5115 A) and the red (6300 — 7250 A) 
bands, giving a resolution of 1.3 A FWHM. We list our observations in Table 2. 
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3. PRE-OUTBURST STATE (1996 APRIL) 
3.1. Overview 

GRO J1655— 40 was in quiescence between late 1995 August and 1996 April. In 1996 March, 
its soft X-ray luminosity, inferred from the ASCA observations (at 2 — 10 keV), was ~ 2 x 10'^^ erg 
s~^ (see Orosz &: Bailyn 1997). The source was not detected by RXTE/ASM (at 2 — 12 keV) above 
the intensity level of ~ 12 mCrab before 1996 April 25 (Levine et al. 1996; Remillard et al. 1996); 
it was also undetected by BATSE before May. On April 25.38 ±0.78 UT (HJD 2450198.88 ±0.78), 
the soft (2 — 12 keV) X-ray intensity began to rise (Remillard et al. 1996). The optical brightening, 
however, had already started on April 20 (Orosz et al. 1997). The fitted times of the initial rise 
were April 19.25 ± 0.29 UT (HJD 2450192.76 ± 0.29) for the / band, April 19.37 ± 0.26 UT (HJD 
2450192.88 ± 0.26) for the R band, April 19.82 ± 0.15 UT (HJD 2450193.32 ± 0.15) for the V band 
and April 20.34 ± 0.18 UT (HJD 2450193.84 ± 0.18) for the B band. 

Our spectroscopic observations were conducted between HJD 2450194.26 and HJD 2450194.30 
(April 20), and between HJD 2450195.10 and HJD 2450195.15 (April 21), just after the initial rise 
in the optical brightness, but before the RXTE/ASM detected the rise in the X-ray intensity. On 
the first night we took a series of five 600-s spectra and one 1200-s spectrum (although focussing 
problems on the red arm of the spectrograph affected the quality our red spectra). The binary 
phase during the observations was 0.485 < < 0.501 (ephemeris of Orosz & Bailyn 1997), and 
the inferior conjunction of the secondary star was at 4> = 0.75. On the following night, wc took 
seven 600-s spectra, at the binary phase 0.804 < cj) < 0.824. (If we adopt the ephemeris of van der 
Hooft et al. (1998), the observations were made at phase 0.480 < (f) < 0.496 on the first night, and 
0.799 < < 0.819 on the second.) 

3.2. General spectral features 

The H I Balmer and Paschen lines are clearly seen in absorption in our spectra (Figures 2 and 
3). The He II A4686 emission line, the Bowen fluorescence N III A4634 line and the blend N III 
AA4641,4642 were not detected on either night. These emission lines, which are often seen from 
X-ray irradiated accretion disks (Smak 1981; Horne k. Marsh 1986), were, however, observed 20 
days later (on May 11) by Hynes et al. (1998b). 

By comparing our spectra with the spectra of F giant and subgiant stars and with the spectra of 
the system taken by Orosz & Bailyn (1997) during a quiescent state in February 1996, we conclude 
that the optical continuum emission was predominantly from the secondary star. This conclusion is 
supported by the agreement between the observed velocity shifts of the Balmer lines in our spectra 
and the radial velocity of the secondary star expected at those phases (not shown). 
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4. HIGH-SOFT STATE (1996-1997) 

4.1. 1996 June observations 

In late 1996 April, a transition from quiescence to a high/soft X-ray state occurred. The 
soft (2 — 12 keV) X-ray flux (measured by RXTE/ASM) increased rapidly, reaching ^ 1.5 Crab 
and remaining above that level for more than six months. Around 1996 November the soft X-ray 
intensity began to decline; but in 1997 January it increased again (Figure 1). It then stayed at the 
1.3 Crab level till 1997 July, after which the system entered another quiescent state. 

At the initial stage of the outburst the X-ray spectrum consisted of a blackbody component 
and a truncated power-law tail. The black body temperature was T ~ 0.8 keV. The photon index 
of the power-law tail was F ~ 3.5 and the exponential cutoff was at 2± 13 keV (Hynes et al. 1998b). 
Around 1996 May 27 - 28 (HJD 2450230 - 31) the X-ray spectrum hardened, and the source 
became detectable by BATSE (Figure 1). The photon index of the power-law was F ~ 2.5 on June 
20, and there was no obvious sign of a high-energy cutoff below ~ 100 keV (Hynes et al. 1998b). 
A radio flare was observed by the Molongo Observatory Synthesis Telescope (MOST), apparently 
coincident the X-ray spectral transition (Hunstead &; Campbell- Wilson 1996; Hunstead, Wu & 
Campbell- Wilson 1997). 

The optical continuum had increased a few days before the onset of the soft X-ray outburst 
(see §3.1). After reaching a peak in mid-May, the optical brightness then declined steadily, with 
an e-folding time of about 50 days (Hynes et al. 1998b). By August 1996 it had returned to the 
previous quiescent level. Flickerings on time-scales of a few seconds, apparently correlated with 
X-ray variability, were reported by Hynes et al. (1998a). 

Our observations were carried out on June 8-12 and June 17, during the soft X-ray outburst 
(Figure 1). We obtained fifty simultaneous red and blue spectra, each of 2000 s duration. In 
Figure 4 we show two blue and red spectra taken on June 8 and June 10. The intensity of the 
optical continuum in May was slightly higher than in June (cf. Hynes et al. 1998b). The soft X-ray 
flux was similar in these two epochs; however, the hard X-rays turned on after May 27. 

4.2. General spectral features in the 1996 June observations 

The Ha and H/3 lines showed a broad absorption trough in 1996 June. The absorption com- 
ponent was stronger on June 8 (EW ~ -6 A, FWHM ~ 70 A ~ 3000 km s'^ for Ha; EW ~ -10 
A, FWHM ~ 70 A ~ 4000 km s~^ for H/?) than on the other nights. The absorption troughs 
were partly filled by narrower emission components (Figure 4), which were not seen in the spectra 
obtained by Hynes et al. (1998b) four weeks earlier. 

The equivalent width of the Ha emission seemed to increase with the hard X-ray flux (Figure 5). 
The emergence of the Ha emission component after the May 27 turn-on of the hard X-ray flux also 
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suggests a direct association between these emissions. We have calculated the discrete correlation 
function, defined by Edelson &: Krolik (1988), for the EW of the Ha emission line and the hard 
X-ray flux and have found a 2-cj peak in the correlation function at the lag of 0.0 it 0.5 d, for a 0.5 
d interval bin (Figure 6). A similar correlation is also found in our 1994 August - September data 
(see §6.4.1). 

The average Ha emission profile was double-peaked, with FWHM ~ 900 km s~^ , and peak- 
to-peak velocity separation = 500 ± 50 km . However, the relative strength of the two peaks 
varied substantially over each orbital cycle (see §4.3). 

The He II A4686 emission line had a double-peaked profile (Figure 7) with a peak-to-peak 
separation of 8.5 it 0.5 A, corresponding to a velocity separation of 545 it 30 km . A double- 
peaked He II A4686 emission line had been detected before the X-ray hardening (see Table 6 in 
Hynes et al. 1998b): archival spectra from the Anglo- Australian Telescope show that the peak-to- 
peak velocity separation was 480 it 40 km s~-^ on 1996 May 11. The central position of the line 
varied sinusoidally, with the variations consistent with the projected radial velocity of the primary 
star for a mass ratio q = 0.33 (Soria et al. 1998), suggesting that the emission originated from the 
accretion disk. The average line center was blue-shifted by 40 it 10 it 5 km s"^ with respect to the 
systemic velocity determined by Orosz & Bailyn (1997), where the first source of error is in the 
determination of the line center, the second is the systematic error in the wavelength calibration. 
We have not found any significant correlation between the EW of the He II A4686 line and the hard 
X-ray flux (Figure 5). 

The He I A6678 emission line appeared to be double-peaked, with only one peak visible at 
some orbital phases. The velocity separation, which was determined when the two peaks were 
present, varied between 520 and 650 km over the time of our observations. These velocities are 
consistent with those of the Ha and the He II A4686 lines. The EW of the line was about 0.7 A. 

A weak narrow absorption component was also detected at A6678 (not shown), superimposed 
on the broader, double-peaked emission component. The narrow absorption component was visible 
at most orbital phases and had an EW ~ —0.2 A and an FWHM 2± 4 A (2± 170 km ). A similar 
narrow absorption component was also found at A7065. 

The velocity shifts of the narrow absorption components of the He I lines were consistent with 
the projected radial velocity of the secondary star (Figure 8), suggesting that these lines were 
due to absorption near the secondary star. The apparent small discrepancy between the observed 
and predicted velocities may be due in part to the uncertainty in our measurement of the line 
position, in part to additional absorption by cold gas near the rim of the accretion disk or in the 
accretion stream, or to a non-uniform absorption by the photosphere of the secondary star, which 
was strongly irradiated on the side facing the primary. 

We note that He I A7065 showed a stronger broad absorption component briefly on June 8, 
with EW ~ —1.8 A and FWHM ~ 35 A ~ 1500 km s"^ . It disappeared later on the same night 
while the hard X-ray flux increased, and was not seen during the following nights (see Figure 4). 
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Narrow, single-peaked N III A4634 and N III AA4641, 4642 Bowen fluorescence lines were ob- 
served in 1996 June. As these lines were also seen before the hard X-ray turned on [see the spectra 
taken by Hynes et al. (1998b) in 1996 May], they were probably uncorrelated with the hard X-ray 
flux. By comparing the Bowen line profiles and radial velocity shifts with those of the He II A4686 
line, we conclude that the Bowen fluorescence lines and the He II line originated from different 
regions. Judging from the radial velocities shifts of the Bowen lines, Soria et al. (1998) suggest 
that they were emitted from a localized region in the outer accretion disk, possibly a hot spot in 
phase with the secondary star. 

4.3. Profile of the Ha emission component 

Among the emission lines observed in 1996 June, the Ha line had the highest signal-to-noise 
ratio. If we adopt a reddening E(B — V) = 1.2 (Hynes et al. 1998b), the inferred average flux for 
the Ha emission line was ~ 5 x 10~^^ erg cm~^ s~^ (cf. -Fue ii ~ 3 x 10~^^ erg cm~^ s^^ for 
the He II A4686 line). Although the EW of the Ha line appeared to be correlated with the hard 
X-ray flux (Figure 5), we have found that the velocity proflle of the line, with its peak normalized 
to unity, depended only on the binary phase. 

The systematic variation of the Ha emission line profiles over an orbital cycle can be seen 
in the sequence of profiles plotted in Figures 9, 10 and 11. The red-shifted peak was stronger 
than the blue-shifted peak between phase 0.25 and 0.85, with the blue-shifted peak dominating at 
other phases. By considering the averaged, normalized Ha line profiles over the whole epoch of our 
observations (see Figure 2 in Soria et al. 1998), and the averaged profiles obtained between orbital 
phases 0.80 and 0.92 and between phases 0.19 and 0.28 (Figure 12), when the highest signal-to-noise 
ratio was achieved and both peaks had comparable strength, we deduce that the FWHM of the line 
was ~ 900 km s^^ , and the peak-to-pcak velocity separation was 500 it 50 km s~^ . These values 
are similar to those determined for the double-peaked He I and He II emission lines. 

In principle, sjnnmetric double-peaked lines can be produced by an accretion disk (Smak 1981; 
Horne & Marsh 1986). As the Ha line profile changed with a period equal to the binary period, 
this eliminates the possibility that the line asymmetry was due to a precession of the accretion 
disk, which should have a period longer than the binary period (e.g. Kumar 1986, Warner 1995). 
Neither could it be caused by the secondary star eclipsing either side of the disk, because our data 
show asymmetry of the two peaks even when the star was behind the accretion disk. 

If the emission were from the irradiatively-heated surface of the secondary star, we would 
expect a stronger blue component in the emission line when the star was approaching (at phases 
0.25 — 0.75), and vice versa. This is inconsistent with the observation that the blue peak of the 
line was stronger at around phase 0. Therefore we reject the possibility that the secondary star was 
the dominant source of the emission. (See §5 for more detailed discussion on the Ha emission.) 
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4.4. 1997 June observations 

The soft X-ray flux in 1997 June was slightly lower than in 1996 June, while the hard X-ray 
flux was signiflcantly lower (Figure 1). We observed the system on 1997 June 14 and 15. No optical 
photometric observations were carried out during this period. Although the conditions were not 
photometric, a comparison with other stars in the fleld led to an estimate ofV ^ 16.5 for the source, 
about half a magnitude fainter than in 1996 June. 

The Ha emission line had a red-shifted peak around phase 0.5 and a blue-shifted peak around 

phase 0, as in 1996 June; the signal-to-noise ratio was much lower than in 1996. The velocity shifts 
with respect to the center of mass of the binary were ~ 255 km for the red-shifted peak, and 
~ —230 km s^^ for the blue-shifted peak, similar to those measured in 1996. 

For the Ha hne, EW = 2.5 ± 0.4 A on June 14, and 2.7 ± 0.3 A on June 15. The emission 
component of the H/3 and higher Balmer lines was not detectable. The Balmer emission was 
therefore weaker than in 1996 June. The He II A4686 line was seen in emission with EW = 3 ± 1 
A, as strong as in 1996. The Bowen fluorescence N III lines were also detected. 

5. BALMER EMISSION IN THE HIGH-SOFT STATE 
5.1. Kinematic interpretation of the Ha emission 

5.1.1. Peak separation 

Figure 13 shows the velocity shift of the peaks of the Ha emission component detected in 
the 1996 June and the 1997 June spectra. Only one peak was present at some orbital phases (cf. 
Figures 9, 10 and 11). The velocity separation of the peaks (~ 500 km s~^ ) was larger than the 
radial velocity amplitude of the secondary star. This implies that the emission region was located 
within the orbit of the secondary star, thus ruling out the possibility of emission from the surface 
of the secondary (see also §4.3) or from a circumbinary disk. 

According to Orosz &; Bailyn (1997), the projected radial velocity semi-amplitude of the sec- 
ondary star in GRO J1655— 40 is K2 = 228.2ib2.2 km s~^ , and the masses of the two components are 
Ml ~ 7Mq and M2 ^ 2.3Mq. If we assume a thin, non-circular, Keplerian accretion disk truncated 
at the tidal radius, as described in Paczyhski (1977), then the disk emission lines would have peaks 
with observed radial velocities at each phase given by the dash-dotted lines in Figure 13. The tidal 
truncation radius is approximately given by = 0.60 a/(l -|- q) for 0.03 < q < 1 (Warner 1995), 
where a is the separation between the centers of mass of the binary components. The peak-to-peak 
velocity separation averaged over an orbital cycle would be ~ 820 km . Phillips et al. (1999) 
derived lower values for the radial velocity semi-amplitude of the secondary star and for the masses 
of the binary components. They obtained 192 < K2 < 214 km , 4.1 < Mi < 6.6M0 and 
1.4 < M2 < 2.2Mq (90 per cent confldence limit). For a fixed mass ratio q = 1/3, this implies 
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that the peak-to-peak velocity separation averaged over an orbital cycle, expected from a tidally- 

truncated Keplerian disk, would be 690 < AV < 760 km . Peak separations ~ 770 km s~^ would 
instead be expected if we adopt the parameters determined by Shahbaz et al. (1999), who inferred 
a velocity semi-amplitude K2 = 215.5 it 2.4 km s^^ for the companion star. 

The velocity separations deduced from our 1996 June and 1997 June data are lower than the 
predicted values, especially for the red-shifted peak at phases around 0.5 and the blue-shifted peak 
at phases around 0. Therefore, the double-peaked profile of the Ha emission component cannot be 
explained by a conventional thin Keplerian accretion disk truncated at the tidal radius. We note 
that the peak-to-peak velocity separations of He II A4686 and He I A6678 were also too small to be 
consistent with a tidally-truncated Keplerian disk. 

Remarkably low rotational velocities, inconsistent with Keplerian disks truncated at their 

tidal radii, have also been inferred from the double-peaked Ha emission line profiles observed from 
A0620— 00 and GS 1124—68 (Orosz ct al. 1994). A possible explanation for the low rotational 
velocities inferred for the outer edge of the accretion disks in these three systems is that the disks 
extended slightly beyond their tidal radii, (cf. Whitehurst 1988). 

5.1.2. A schematic model 

We have found that the observed variations in the Ha line profile can be explained by the 
model shown in Figure 14, with two separate emission regions on the accretion disk. The model 
can reproduce the low radial velocities of the peaks and their behavior at various orbital phases. 
For example, only the red-shifted component would be visible around phase 0.5 and only the blue- 
shifted component around phase 1, while both peaks would be visible around phases 0.25 and 
0.75. 

Non-axisymmetric patterns of emission similar to our model were proposed by Steeghs, Har- 
laftis (fe Horne (1997) to explain the Balmer and helium emission line profiles observed from the 
dwarf nova IP Peg, and by Neustroev & Borisov (1998) to explain the Balmer line profiles observed 
from the dwarf nova U Gem. In both cases, they were interpreted as evidence of a two-armed spiral 
density wave or shock in the accretion disk, induced by the tidal interaction with the companion 
star. In the case of IP Peg, the observed rotational velocities also suggested that the disk could 
extend beyond its tidal radius (Steeghs et al. 1997). 

It has been suggested that accretion stream overflows, hot spots, and uneven illumination of the 
accretion disk from the central X-ray source can also produce an anisotropic emissivity. However, 
it is beyond the scope of the present paper to explore in details why the Balmer emission regions 
might assume the geometrical configuration shown in our model. 

Furthermore, a narrow Ha absorption component from the companion star is also likely to be 
present, as suggested by the detection of other stellar absorption features from He I and Fe I. It 
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can be noticed from the radial velocity curve of the companion star that the presence of a stellar 
absorption line superimposed on the disk emission would also be qualitatively consistent with the 
observed behavior of the line profile at various orbital phases, if absorption and emission have 
comparable strength (Figures 9, 10 and 11). However, we argue that the effect of anisotropic disk 
emission is more significant than the effect of a stellar absorption line, and can better explain the 
velocity shifts of the emission peaks observed over an orbital phase. 

One might expect the Ha emission line profiles to be roughly symmetric at superior and 
inferior conjunction; however, the Ha line in general showed a red-shifted peak stronger than the 
blue-shifted peak, and a blue wing more extended than the red wing (see Figure 9 and Figure 11). 
The asymmetry of the profiles may indicate the presence of both opacity and kinematic effects, 
and may be evidence of a thin disk-wind. It is worth noting that P-Cygni profiles were observed 
in some UV resonance lines from the system (Hynes et al. 1998b). 

5.2. Formation of absorption and emission lines 

5.2.1. Broad absorption lines 

The Ha and H/3 lines in the spectra obtained during the 1996/1997 outburst all show broad 
(FWHM ~ 3000 km s^^ ), shallow absorption components. Broad absorption components are also 
seen in the 1994 August - September spectra, when the system was in outburst (see §6.2), but not 
in the 1996 April spectra, when the system was in quiescence. 

Absorption lines can be produced in an accretion disk if the disk is optically thick and its 
temperature decreases with height above the central plane. The absorption features seen in our 
spectra are broader than the emission cores: this suggests that they were probably formed in the 
inner part of an optically thick disk, where viscous heating from the central plane dominates over 
external irradiative heating. 

Broad, shallow absorption features at the Balmer lines have been seen in the spectra of other 
BHCs in outburst, such as A0620-00 (Whelan et al. 1977), GS 1124-68 (Delia Valle, Masetti & 
Bianchini 1998) and, most prominently, GRO J0422+32 (Casares et al. 1995). Similar features are 
also often present in the spectra of UX UMa stars (e.g. Warner 1995) and in the outburst spectra 
of Dwarf Novae (Robinson, Marsh &; Smak 1993). 

5.2.2. Double-peaked emission lines 

Double-peaked lines can be emitted from a hot temperature-inversion layer on the X-ray irra- 
diated surface of the accretion disk. Provided that the spectrum is sufficiently soft. X-rays can be 
absorbed at a small depth, forming a thin temperature-inversion layer but leaving most of the ver- 
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tical structure of the disk undisturbed (see Tuchman, Mineshige & Wheeler 1990; Wu et al. 1999). 

The strong, double-peaked He II A4686 line seen during the high-soft state in 1996-1997 was 
probably emitted via radiative recombination in this thin, hot layer, at temperatures ~ 10^ K. 

While the He II A4686 was detected throughout the 1996 outburst, Balmer emission was seen 
only after the hard X-rays turned on. Balmer lines are likely to be emitted at lower temperatures 
{•-^ 10^ K). Moreover, the Ha double-peaked profile is strongly asymmetric and phase-dependent 
(§5.1.2), while the He II A4686 profile is approximately symmetric at all phases. For all these 
reasons we suggest that He II A4686 and the Balmer lines were emitted from different regions: the 
Balmer lines probably originated from a deeper, denser layer in the disk where matter was heated 
by harder X-ray photons. 

6. ONSET OF A HARD X-RAY OUTBURST (1994 AUGUST - SEPTEMBER) 

6.1. Overview 

GRO J1655— 40 was active in the radio and hard X-ray energy bands in 1994 August - Septem- 
ber. The 843 MHz radio flux density, measured by the MOST, was declining after a large outburst 
that reached about 8 Jy in early August (Wu &; Hunstead 1997). On 14 September another radio 
outburst started. It was weaker than the previous outburst, with a peak flux density of about 2 Jy. 

The hard X-ray flux recorded by BATSE (Harmon et al. 1995) was below the 0.7 photon 
cm~^ s"^ level during the period August 15 - September 4. A sharp rise in the X-ray flux occurred 
on September 5, and the flux stayed at the 2 photon cm~^ s~^ level for about 10 days. It then 
declined abruptly, apparently in coincidence with the rise in the 843 MHz radio luminosity around 
September 14 (see Fig. 1 in Wu &; Hunstead 1997). 

The span of our observations covered the transition phase around the hard X-ray rapid increase. 
Good signal-to- noise spectra, with a resolution of 1.3 A FWHM, in the spectral regions centered 
at Ha and at He II/ H/3, were obtained on August 30 - September 4, before the rise in the hard 
X-ray flux. The sampling of binary phase was roughly uniform over the six- night observations. Two 
lower-resolution spectra in the Ha and He II/ H/? regions were taken on 1994 September 6, after 
the rise in the X-ray flux, and showed a dramatic change compared with the August 30-September 
4 spectra. 

6.2. General spectral features before September 6 

In Figure 15 we show the red spectrum, centered at Ha, averaged over the six nights from 
August 30 to September 4. As in the spectra obtained during the 1996 - 1997 high-soft state, we 
notice a broad absorption component partly filled by narrow Ha emission. The average EW of the 
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absorption component is -6 ± 1 A, and its average FWHM 60 ± 10 A (~ 2700 ± 500 km s~^ ). The 
narrow emission component has an average EW = 5.3 it 0.2 A and an average FWHM = 450 it 20 
km s~^ . 

Figure 16 shows the averaged blue spectrum, centered at He II/ H/3. The broad absorption 
component of H/3 has an EW = -4 ± 1 A and an FWHM = 50 ± 10 A (~ 3000 ± 600 km s^^ ). A 
narrow emission component is superimposed on the broad absorption trough, with EW = 1.1 it 0.1 
A and FWHM = 550 ± 30 km 8"^ 

He II A4686 appears in the spectrum as a strong narrow emission line (Figure 16). The EW 
of the line is 5.2 ± 0.2 A and its FWHM is 540 ± 20 km s"^ . Weaker emission is detected at He II 
A4542, with EW = 0.5 ± 0.2 A and FWHM = 800 ± 50 km s-^ . 

A broad emission component (FWHM = 950 it 100 km s~^ ) was detected on each night at 
He I A6678. An additional narrower component (FWHM = 450 it 50 km s~^ ) was detected on some 
nights; it was particularly strong on August 30 (cf. Figures 17). We note that this line might be 
contaminated by He II A6683 emission. 

The spectra show a broad, fiat-topped line at about 5005 A. It cannot be attributed to [O III] 

A5007 both because the wavelengths are discrepant, and because we do not see any emission from 
[O III] A4959, whose strength should be one-third of that of [O III] A5007. (The only forbidden 
lines possibly detected in our spectra obtained between 1994 August 30 and September 4 are [N II] 
A6548, blended with the much stronger blue wing of Ha, and [N II] A6589; see Figures 17). We 
identify the line at about 5005 A as N II A5005. 

Another strong low-ionization metal line detected in the spectrum is the blend O II AA4941, 4943. 
Other broad emission lines detected from N II and O II are listed in Table 3. We do not include a 
weak emission line at 6375 A in Table 3, as it was seen only on the first two nights (EW = 0.7ib0.1); 
its wavelength is consistent with that of an Fe X line. Bowen N III emission lines were also observed, 
broader than in 1996 June, and with a radial velocity curve consistent with that of the primary 
rather than the secondary. 

6.3. Line classification 

The FWHMs of the optical lines observed before 1994 September 6 span a large range of 
values, from ^ 400 km s~^ to ~ 3000 km s~^ . The broadest lines are those seen in absorption at 
Ha and H/3. We notice that a broad absorption trough at H/3 was also detected in the 1995 March 
outburst by Bianchini et al. (1997), and during the 1996 May-June high-soft state. The physical 
interpretation of the broad Balmer absorption lines is probably the same as discussed in §5.2.1. 
Henceforth, we will focus only on the emission lines. 

For the orbital parameters of GRO J1655— 40, the projected rotational velocity of the outer 
rim of a thin, Keplerian accretion disk, truncated at or slightly beyond its tidal radius, is ^ 350 
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km s~ . Therefore, any lines with an FWHM < 700 km s~ cannot come from a thin Keplerian 
disk. We classify the emission lines observed before 1994 September 6 as "broad" and "narrow" 
(Table 3), according to whether their FWHM was larger or smaller than the minimum value of 
FWHM expected for disk emission. 

The broad emission lines were usually flat-topped, while the narrow emission lines were either 
single-peaked or had a hint of a double-peaked profile but with very low velocity separation. In 
some cases the distinction is blurred, especially for the weak lines. However, in most situations the 
two kinds of line profiles can be discerned. We notice that both the EW and the FWHM of the 
narrow lines changed significantly from night to night, while EW and FWHM of the broad lines 
were more stable. 

Our phenomenological classification of broad and narrow emission lines in GRO J1655— 40 is 
similar to the classification of emission lines in AGNs. For example, a study of 123 high-luminosity 
AGNs (Wills ct al. 1993) has shown that the profile of the C IV A1549 emission line consists of a 
"broad" base and a "narrow" core, which are emitted from distinct regions. The FWHM of the 
line is determined by the core/base ratio: when the core is dominant, the line will appear narrow 
and sharply peaked; and when the base is dominant the line will be broad and flat-topped (cf. Fig. 
2 in Wills et al. 1993). 

6.3.1. Origin of the "broad" emission lines 

The low-ionization metal lines from N II and O II detected in our spectra are examples of 
broad emission lines; their FWHMs are consistent with a disk origin. Unlike the disk emission lines 
seen in 1996 June, their profiles were generally flat-topped rather than double-peaked. 

Flat-topped lines can be emitted in a wind from the surface of the accretion disk. Murray 
& Chiang (1997) showed that if the lines have substantial optical depth (r/ ^ 1), the contribution 
of the wind from the front and back sectors of the disk (where the projected radial velocities are 
lower) is enhanced relative to the contribution of the wind from the sides of the disk (where the 
projected radial velocities are higher), thus smearing out the two peaks in the line profile. 

From the fact that we detected only low-ionization broad lines from He 1, O II, N II and 
N III, we infer an ionization parameter U = Ly^/n^r'^ ~ 10 in the line-emitting disk wind. The 
corresponding temperature would be ~ 1 - 2 x 10^ K if the gas was collisionally ionized (Hatchett, 
Buff k, McCray 1976); the temperature would be lower if the line-emitting gas was photoionized. 

The possibility of an outflow or disk wind suggested by the broad optical emission line profiles 
observed in the 1994 X-ray outburst was already discussed in Hunstcad ct al. (1997), who compared 
the spectrum of GRO J1655— 40 with the spectra of WN6-8 Wolf-Rayct stars. In spite of the 
similarities, there are several differences between the two cases. Firstly, high-ionization emission 
lines usually found in WR stars, e.g. from N IV and N V, were not observed, indicating a lower 
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temperature for the outflow in GRO J1655— 40. Secondly, no forbidden [O III] lines were observed, 
implying a higher density. Thirdly, the broad, flat-topped proflles are better explained by a disk- 
wind model rather than by a spherical outflow as in typical Wolf-Rayet stars. 

6.3.2. Kinematics of the "narrow" emission lines 

The narrow emission lines showed large night-to-night changes in their EWs and line profiles. 
In particular, they showed transitions between narrow, single-peaked and slightly broader, double- 
peaked profiles over the time-span of our observations. We notice that He II A4686 line was generally 
broader and more clearly double-peaked than Ha (Figure 18; cf. also the double-peaked profile in 
the combined spectrum shown in Figure 16); its peak-to-peak velocity separation was however 
always < 400 km s~^ . 

Some lines (most notably He I A6678) were probably emitted from both the broad- and the 
narrow-line regions, with the intensity ratio determined by the relative contribution of the two 
components, i.e. as in the "base" vs "core" classification scheme suggested by Wills et al. (1993) 
in their study of AGNs. 

Both the Ha and He II A4686 emission lines were, on average, asymmetric; they were slightly 
blue-shifted with respect to the central line wavelength (taking into account the systemic velocity). 
In the combined spectrum of all six nights, the Ha line is blue-shifted by 55 ± 9 km s~^ , the H/3 
line by 61 it 12 km s~^ , the narrow component of the He I A6678 line by 82 it 22 km , and 
the He II A4686 line by 41 it 13 km . We note that a similar phenomenon is often observed in 
the emission lines from quasars and in some cataclysmic variables, and can be attributed to the 
expansion of clouds or to a gas outflow, when the receding gas is occulted from our view. 

6.3.3. Optically thin narrow-emission-line region 

Owing to the well-determined orbital parameters, the upper limit to the size of the accretion 
disk in GRO J1655— 40 is reasonably well constrained. We have found that the velocities inferred 
from the narrow widths of the Ha, H/? and He II A4686 lines were signiflcantly lower than the 
velocities expected from a Keplerian disk contained in the Roche lobe of the compact star. 

One might argue that the low FWHM of the narrow emission lines was due to their formation 
in a circumbinary disk (Artymowicz & Lubow 1996). If this were true, no significant velocity 
shifts from the center-of-mass systemic velocity should be seen over an orbital cycle. However, 
Soria et al. (1998) have measured the radial velocity shifts of the line wings at one-quarter of 
the maximum intensity above the continuum, and found a sinusoidal variation consistent with the 
one observed in 1996 June, which was interpreted as the modulation due to orbital motion of the 
compact star. Moreover, the profiles of these lines changed substantially from night to night. In 
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the case of Ha, the line profile was single-peaked on four nights, but double-peaked on the other 
two, with a velocity separation larger than expected from a circumbinary disk. We can therefore 
rule out the possibility of emission from a circumbinary disk. 

We interpret the narrow lines, instead, as emission from a spheroidal or extended optically 
thin region high above the accretion disk (at lower rotational velocity). We notice that some 
cataclysmic variables show narrow, often single-peaked Balmer and He II emission lines, yielding 
rotational velocities much slower than Keplerian. For example, there is strong evidence of an 
optically thin line-emitting region extended well above the disk plane in the eclipsing system PG 
1012-029 (Honeycutt, Schlegel h Kaitchuck 1986). The larger FWHM of the high-ionization He II 
line suggests that it was probably emitted closer to the disk plane than Ha, and/or at slightly 
smaller radii, where we expect the irradiation from the central source to be stronger. 

The night-to-night variability in the line profiles (Figure 18) might reflect changes in the 
line-of-sight opacity of the optically thin region. Opacity variations can be caused by matter 
inhoniogeneities (i.e. moving clouds) in the line emission region or variable obscuration along the 
light path. Furthermore, the variability between single-peaked and double-peaked profiles indicates 
that the extended line-emitting region was not in steady state; for instance, the profiles would 
tend to become broader and more "disk-like" (as those observed in 1996 June) as the thin clouds 
dissipated, and less "disk-like" when the emission from the clouds was stronger. 

The presence of optically thin line-emitting gas well above the disk plane may be explained 
in various ways: some accreting matter might not yet have settled down onto the disk, after the 
end of the strong hard X-ray flare observed by BATSE in 1994 August. We shall argue in §6.4.2 
that the outer disk is likely to be disrupted and puffed up into a thick torus during a hard X-ray 
flare. Some of the gas in the extended narrow-line region may also have been the result of bipolar 
outflows from the disk surface, as suggested by the systematic blue-shift and by the detection of 
strong winds in the broad-line region. 

The extended, optically thin narrow-line-region model is valid only if the UV and soft X-ray 
flux is not too high. Soft irradiating flux of order of the Eddington limit at distances R < 10^^ cm 
would yield an ionization parameter U > 10^ in the thin narrow-line region, for a gas with densities 
nu < 10^^ cm^^. This would imply temperatures ^ 10^ K, difficult to reconcile with the presence 
of lines such as Ha. Although no UV and soft X-ray data were obtained at that epoch, the fact 
that we detected only low-ionization metal lines (for example from N II and O II) implies that the 
soft flux could not have been too strong in the 1994 X-ray outburst. 
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6.4. Spectral transition on 1994 September 6 

6.4- 1- General features 

A sharp increase in the hard X-ray flux was detected by BATSE around September 5 (Harmon 
et al. 1995). The broad N II, N III and O II emission hnes disappeared from the blue spectrum 
obtained on September 6 (Figure 19), and He II A4686 and the Bowen N III lines were much weaker 
than in the previous nights. Ha and H/? were instead much stronger: for Ha, EW = 100 it 5 
A, and for H/?, EW = 8.0 ± 0.5 A. The low-resolution of the red spectrum does not allow for a 
reliable determination of the FWHM of Ha; for H/3, FWHM = 530 ± 50 km s"^ . The Ha/H/? EW 
ratio increased by a factor of 4 between September 4 and September 6. Other lines clearly seen in 
emission are the Paschen series of H I and several He I lines (A4922, A5016, A5876, A6678, A7065). 
All these emission lines were single-peaked; unfortunately, the low S/N of the blue spectrum and 
the low resolution of the red spectrum do not permit an accurate determination of their strengths 
and widths. 

The sudden increase in the Balmer emission between September 4 and September 6 coincided 
with the sharp increase in the hard X-ray flux. The EWs (in absolute value) of Ha, 11(3 and 
He II A4686 are shown in Figure 20 compared with the hard X-ray flux measured by BATSE. The 
Balmer EWs and the hard X-ray flux appear well correlated, as in the 1996 June observations. 
Such correlation is not seen for He II A4686. In fact, its intensity dropped when the the hard X-ray 
flux rose. 

Shrader et al. (1996) also noticed the dramatic increase in the Ha emission and the disappear- 
ance of He II A4686 around TJD 9600, after the hard X-ray flux surge (see their Fig. 1). Tlicir 
low-resolution spectrum obtained on September 8 showed only emission from the low-ionization 
lines H I, He I and O I. 

Optical spectra taken by Bianchini et al. (1997) on 1994, August 12 - 15, during the decline of 
an earlier hard X-ray flare, showed emission almost as strong as that in our September 6 spectrum 
from low-ionization lines like Ha (mean EW = 68.4 A) and H/? (mean EW = 7.4 A). The Ha/H/3 
EW ratio was 9.2, a factor of 3 higher than on September 4, and slightly lower than on September 
6. The high-ionization He II A4686 (mean EW = 5.3 A) and N III AA4641, 4642 (mean EW = 5.9 A) 
lines were also present, showing disk or disk-wind signatures (see Fig. 2 in Bianchini et al. 1997). 
From the peak separation and the width of the flat-topped profiles, we can infer that the outer-disk 
radius was ^ (1.3 ± 0.1) x 10^^ cm. Broad absorption troughs were not seen in those spectra. 

We suspect that the conditions on August 12 - 15 were intermediate between those we observed 
before and after September 5. In the period of low X-ray activity between August 15 and September 
5, a geometrically thin accretion disk was probably re-forming and re-expanding to its tidal radius 
after having been disrupted in the previous hard X-ray flare, only to be then "evaporated" again 
into an extended atmosphere or cocoon by the following flare. 
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6.4-2. Optically thick Balmer emission 

As the photoionization cross-section for hydrogenic ions falls approximately as v^'^l'^ for photon 
energies larger than the threshold energy (13.6 eV for hydrogen) (see Osterbrock 1989), the hard 
(20 — 100 keV) X-ray photons contribute little to the ionization level of the gas, compared with 
the UV and soft X-ray photons. However, Compton heating by the hard X-rays can alter the 
disk structure and drive a disk-wind (e.g. Begclman, McKee & Shields 1983; Man & Shaviv 1996). 
In the extreme conditions, the dense disk-wind would create a cocoon surrounding the accretion 
disk; in this situation one may consider the accretion disk being practically evaporated into a thick 
atmosphere. The absence of the characteristic broad absorption troughs at Ha and H/3 and of 
broad disk-like emission lines in our September 6 spectrum suggests the absence of an exposed, 
geometrically thin, optically thick accretion disk. 

An increase in the optical depth of the corona could have caused the increase in Balmer line 
emission and the disappearance of the high-ionization lines. The higher Ha/H/? EW ratio measured 
during the 1994 August and September hard X-ray flares, and the Paschen emission lines observed 
on September 6, indicate that the corona was optically thick in the Balmer series. When this 
situation occurs, every H/3 photon is resonantly absorbed by a nearby hydrogen atom in the n = 2 
state; the process of emission and reabsorption is repeated until the H/3 photon splits into Ha -|- 
Pa photons, which will both escape ("Case C" recombination, see e.g. McCray 1996). 

After the end of the 1994 August hard X-ray flare, the extended cocoon probably became 
optically thin to the Balmer lines ("Case B" recombination), as observed on August 30 - September 
4. The disk was "puffed up" again by the hard X-ray flare after September 4, so that Balmer lines 

were emitted via Case C recombination. Archival spectra obtained from the Anglo-Australian 
Observatory, taken on September 27 (Figure 21), show features almost identical to those observed 
on August 30 - September 4 (cf. Figures 15 and 16), suggesting that the cocoon once again became 
optically thin after the end of the second hard X-ray flare. 

The simultaneous decrease in the observed strength of the He II and Bowen N III lines at the 
time of the hard X-ray increase can perhaps be attributed to a geometric effects. As GRO J1655— 40 
has a high orbital inclination, a thick outer-disk rim, inflated by hard X-ray irradiation, may occult 
the inner region where the high-ionization lines are formed, and explain their weakness in the 
September 6 spectrum. 

Finally, we note that a major episode of plasma ejection occurred on 1994 September 9, inferred 
from the radio observations (Harmon ct al. 1995; see also Wu & Hunstead 1997), a few days after 
the change in the optical and hard X-ray spectra. It was suggested (Fabian & Rccs 1995) that there 
may be a connection between the production of collimated jets and the presence of an extended hot, 
high pressure atmosphere and a quasi-spherical accretion flow in the centers of elliptical galaxies, the 
usual hosts for radio-loud emission. Whether the same mechanism was a catalyst for the collimated 
ejections observed in the radio-loud outburst of GRO J1655— 40 in 1994 August - September is an 
issue worth investigating. 
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7. SUMMARY 

We have obtained optical spectra of the soft X-ray transient GRO J1655— 40 during different 
energy states (quiescence, high-soft and hard outburst) between 1994 August and 1997 June. 

Spectra taken one day after the observed brightening in V in 1996 April did not at that stage 
show evidence of emission lines from an accretion disk. The main spectral features were narrow 
absorption lines from the secondary star. The optical brightening was probably due to an increase 
in the continuum emission from the disk; emission lines could not be formed at that stage, in the 
absence of X-ray irradiation. 

For the 1996 - 97 high-soft state, we have identified characteristic features such as broad 
absorption lines at Ha and H/3, partly filled by double-peaked emission lines. Wc argue that the 
broad absorption was formed in the hot inner disk, and the double-peaked He II A4686 emission 
originated in a temperature-inversion layer on the disk surface, created by the soft X-ray irradiation. 
The observed rotational velocities suggest that the disk was probably extended beyond its tidal 
radius. We have also found that the double-peaked Ha emission was associated with the hard X-ray 
flux, suggesting that it was probably emitted at comparable radii but from deeper layers (at higher 
optical depth) than He II A4686. We note that the Balmer lines were not emitted uniformly from 
the whole disk surface, but appeared to come only from a double-armed region, possibly the effect 
of tidal density waves or shocks on the disk. 

Wc have identified three classes of lines in the spectra taken in 1994 August - September 
before the onset of a hard X-ray flare: broad absorption, broad (flat-topped) emission and narrow 
emission. The narrow emission line proflles cannot be explained by a conventional thin accretion 
disk model. We speculate that the system was in a transient state in which the accretion disk 
had an extended, optically thin cocoon and significant matter outflow, which would also help to 
explain the systematic blue-shift of the narrow emission lines and the flat-topped proflles of the 
broad emission lines. 

After the onset of a hard X-ray flare, the disk signatures disappeared, and strong Ha and 
Paschen emission was detected, suggesting that the cocoon became optically thick to the Balmer 
lines. High-ionization lines disappeared or weakened. Two weeks after the end of the flare, the 
cocoon appeared to be once again optically thin. 
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Table 1. Physical parameters for GRO J1655— 40 



Ref. 



Distance 


d = 3.2± 0.2 kpc 


a 




3 < d < 5 kpc 


b 


Binary period 


P = 2.62157 ±0.00015 d 


c 




P =2.62168(14) d 


d 


Inclination angle 


i = 69°. 5 ±0°. 08 


c 




63°. 7 < i < 70°.7 


d 


Mass ratio 


q = 0.334 ± 0.009 


c 




0.337 <q< 0.436 


e 


Mass function 


fx = 3.24 ± 0.09 Mq 


c 




fx = 2.73 ± 0.09 Mq 


e 




1.93 < fx< 2.67 Mq 


f 


Mass of primary 


Ml = 7.02 ± 0.22 Mq 


c 




5.5 < Ml < 7.9 Mq 


e 




4.1 < Ml < 6.6 Mq 


f 


Mass of secondary 


M2 = 2.34 ±0.12 Mq 


c 




1.7 < M2 < 3.3 Mq 


e 




1.4 < M2 < 2.2 Mq 


f 


Systemic velocity 


7 = -142.4 ± 1.6 km s"^ 


c 




7= -141.9 ±1.3 km s"^ 


e 




-153 < 7 < -143 km 8"^ 


f 


Quiescent V mag 


17.3 


g 



^Hjellming &i Rupen 1995 

''Tingay et al. 1995 

'^Orosz k Bailyn 1997 

'^van der Hooft et al. 1998 
^Shahbaz et al. 1999 
^Phillips et al. 1999 
sBailyn et al. 1995 
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Table 2. Our spectroscopic observations of GRO J1655— 40 



Date 


Epoch 


Phase range 


Wavelength range 


Resolution 




(HJD - 2449000) 




(A) 


(A FWHM) 


AAT 3.9 m 


1994 August 26 


590.969-590.976 


0.359-0.362 


5645-7240 


3.5 


1994 August 29 


593.947-593.954 


0.495-0.498 


3780-7450 


7.3 


1994 August 30 


594.866-594.883 


0.846-0.852 


6280-6825 


1.3 




594.897-594.909 


0.858-0.862 


4432-5051 


1.3 


1994 August 31 


595.863-595.877 


0.226-0.231 


6280-6825 


1.3 




595.885-595.899 


0.234-0.240 


4432-5051 


1.3 


1994 September 1 


596.867-596.883 


0.609-0.615 


6280-6825 


1.3 




596.888-596.911 


0.617-0.626 


4432-5051 


1.3 


1994 September 2 


597.857-597.870 


0.987-0.992 


6280-6825 


1.3 




597.872-597.885 


0.992-0.997 


4432-5051 


1.3 


1994 September 3 


598.860-598.873 


0.369-0.374 


6280-6825 


1.3 




598.875-598.894 


0.375-0.382 


4432-5051 


1.3 


1994 September 4 


599.860-599.873 


0.751-0.756 


6280-6825 


1.3 




599.882-599.895 


0.759-0.764 


4432-5051 


1.3 


1994 September 6 


601.861-601.864 


0.514-0.515 


5500-11000 


20 








3925-5500 


2.5 


ANU 2.3 m 


1996 April 20 


1194.259-1194.303 


0.485-0.501 


3600-9300 


4.8 


1996 April 21 


1195.096-1195.151 


0.804-0.824 


3600-9300 


4.8 


1996 June 8 


1242.980-1243.327 


0.069-0.202 


4360-5320 


1.3 




1242.908-1243.327 


0.042-0.202 


6300-7250 


1.3 


1996 June 9 


1243.893-1244.285 


0.418-0.567 


4150-5115 


1.3 








6300-7250 


1.3 


1996 June 10 


1244.868-1245.292 


0.790-0.951 


4150-5115 


1.3 








6300-7250 


1.3 


1996 June 11 


1245.885-1246.153 


0.178-0.280 


4150-5115 


1.3 








6300-7250 


1.3 


1996 June 12 


1246.877-1247.293 


0.556-0.715 


4150-5115 


1.3 








6300-7250 


1.3 


1996 June 17 


1252.029-1252.052 


0.521-0.530 


6300-7250 


1.3 


1997 June 14 


1614.010-1614.092 


0.599-0.630 


4150-5115 


1.3 








6300-7250 


1.3 


1997 June 15 


1615.056-1615.243 


0.998-0.069 


4150-5115 


1.3 








6300-7250 


1.3 
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Table 3. Main permitted lines detected during the 1994 outburst 



Line 


EW (A) 


FWHM (km s"^ ) 


BROAD ABSORPTION LINES 


H/? A4861 


-4± 1 


3000 ± 600 


Ha A6562 


-6± 1 


2700 ± 500 


BROAD EMISSION LINES 


II A4452 


0.7 ±0.2 


1350 ± 150 


N III A4515 


0.7 ±0.2 


950 ± 50 


He II A4542 


0.5 ±0.2 


800 ± 50 


II AAA4602, 4609, 4613^ 


2.1 ±0.2 


1400 ± 150 


NIII AA464I,4642b 


5.0 ±0.5 


1050 ± 100 


II AA4941,4943'^ 


1.9 ±0.1 


1750 ± 50 


N II A5005 


1.0 ±0.1 


1400 ± 50 


He I Aeers'^'^ 


0.9 ±0.1 


950 ± 100 


NARROW EMISSION LINES 


Ha A6562 


5.3 ±0.2 


450 ± 20 


H/3 A4861 


1.1 ±0.1 


550 ± 30 


He II A4686 


5.2 ±0.2 


540 ± 20 


He I Aeers*^'^ 


0.9 ±0.1 


380 ± 40 



^Possible contribution also from N II AAA4601, 4607, 4614 
and (less likely) from N IV A4606. 

'^after deblending a small contribution from N II A4631; 
possible contribution also from N III A4634. 

'^after deblending a small contribution from He I A4922. 

"broad" component of the line. 

'^Possible contribution also from He II A6683. 

^"narrow" component of the line. 
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Fig. 1. — a: daily averages of tlic liard (20 — 100 keV) X-ray flux measured by BATSE. b: daily 
averages of the soft (2 — 12 keV) X-ray flux measured by RXTE/ASM. The lightcurves cover 
the period from 1996 March 18 to 1997 September 9. The solid lines mark the epochs of our 
observations (1996 April 20-21, 1996 June 8-12, 1997 June 14-15). A transition from quiescence 
to a high-soft X-ray spectral state occurred around MJD 50198 (April 25). The X-ray spectrum 
hardened after MJD 50230 (May 27). A radio flare was observed by MOST on MJD 50232 (May 
28.60 UT) (dashed line). Hjellming (1997) estimates that the peak of the radio flare occurred 
sometime between MJD 50222.5 and MJD 50230.5. The actual onset of the radio flare is not well 
determined: its flux density was already declining with an e- folding time of 1.4 d when the flare 
was detected (Hunstead et al. 1997). In the lower panels, we show an expanded view of the hard 
(c) and soft (d) X-ray flux around the time of our 1996 June observations (shaded regions). 
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Fig. 2. — Flux-calibrated spectra taken on 1996 April 21 (lower spectrum) and on 1996 June 10 
(upper spectrum) with the ANU 2.3 m telescope at Siding Spring Observatory. Wavelengths are 
vacuum heliocentric (cf. Figure 2 in Hynes et al. 1998b). 
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Fig. 3. — Top panel: averaged blue spectrum obtained on 1996 April 20, at spectroscopic phases 
0.485 < (f) < 0.501, according to the cphcmcris of Orosz & Bailyn (1997). Central and bottom 
panels: portions of the averaged red spectrum obtained on 1996 April 21, at spectroscopic phases 
0.804 < (f) < 0.824. Wavelengths are vacuum heliocentric, and the intensities are normalized to 
the continuum. The Balmer series is in absorption. Other strong absorption features are H I 
Paschen series, Ca II AA3934,3968 (the latter is blended with He A3970), and the triplet Ca II 
AAA8542, 8662, 8498 (blended with the Paschen lines). No emission is observed from the N III/C III 
lines, nor from He II A4686. 
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Fig. 4. — a: normalized blue spectrum taken on 1996 June 8, around spectroscopic phase (p = 0.10. 
b: normalized blue spectrum taken on 1996 June 10, around (j) = 0.85. c: normalized red spectrum 
taken on 1996 June 8, around cf) = 0.10. d: normalized red spectrum taken on 1996 June 10, 
around (p = 0.85. Wavelengths are vacuum heliocentric. Broad absorption troughs were detected 
at H/3 and Ha on June 8; they were partly filled on June 10. A broad absorption feature at He I 
A7065 was also present in the June 8 spectrum, but had disappeared on June 10. Emission was 
observed from the Bowen fluorescence lines of N III (single-peaked) and from He II A4686, Ha and 
He I A6678 (double-peaked). The Ha emission was much stronger on June 10 than on June 8. The 
peak-to-peak separations for He II A4686 lines were 525 it 25 km s~^ on June 8, and 545 it 25 km 
on June 10. The peak-to-peak separation for Ha was 550 it 10 km s~^ on June 10; on the same 
night, the peak-to-peak separation for He I A6678 was 650 it 50 km s~^ . 
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Fig. 5. — Top panel: hard (20 — 100 keV) X-ray flux measured by BATSE during the time of our 
1996 June observations. Bottom panel: equivalent width of the Ha and He II A4686 emission lines. 
No correlation was found beween the hard X-ray flux and the EW of He II A4686, but there is 
evidence of a correlation between the hard X-ray flux and the Ha emission. A similar correlation 
was observed in 1994 August - September (cf. Figure 20). 
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Fig. 6. — Discrete cross-correlation function for the equivalent width of the Ha emission line and 
the hard (20 — 100 keV) X-ray flux in 1996 June. The error is large, mainly because of the limited 
number of datapoints available for the correlation function; nonetheless, the two quantities appear 
correlated. Better evidence of a correlation between hard X-ray flux and Balmer emission was 
found in our 1994 observations (cf. Figure 20). 
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Fig. 7. — Normalized line profile of He II A4686, obtained by averaging the spectra taken on 1996 
June 8 at orbital phases 0.08 < cf) < 0.12. Wavelengths are vacuum heliocentric and the intensity 
is normalized to the continuum. The velocity zeropoint is the systemic velocity (7 = —142.4 ±1.6 
km , see Table 1). Weak emission is also detected at positions consistent with C III A4664, O II 
A4699 and O II A4705, as marked in the figure. 
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Fig. 8. — Radial velocity shifts of the narrow absorption components observed at He I A6678 (top 
panel) and He I A7065 (bottom panel) in 1996 June. Dash-dotted line (top panel): best sinusoidal 
fit to the velocity shifts of He I A6678 (semi-amplitude K = 199.9 it 6.4 km s~^ ). Dashed line: 
projected radial velocity of the secondary star according to Shahbaz et al. (1999) (semi-amplitude 
iC = 215.5 ±2.4 km s^^). 
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Fig. 9. — Ha emission line profiles observed on 1996 June 8 (left panel) at spectroscopic phases 
0.06 < < 0.20, and on June 11 (right panel), at phases 0.18 < (p < 0.27. The velocity zeropoint 
(dotted line) is the systemic velocity (7 = —142.4 it 1.6 km ). The projected radial velocities of 
the secondary star with respect to the systemic velocity are also shown next to each spectrum (short 
solid lines). Note that when the secondary star was receding from us, the largest contribution to 
the Ha emission was coming from approaching gas. Both the red and the blue peak were visible 
when the secondary star was near superior conjunction. 
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Fig. 10. — As in Fig. 9, for the 1996 June 9 observations (left panel), at spectroscopic phases 
0.42 < 4> < 0.56, and the June 12 observations (right panel), at phases 0.56 < (p < 0.71. The dashed 
profile in the left panel was obtained on 1996 June 17. Note that while the secondary star was 
approaching us, most of the emission was coming from receding gas. A narrow absorption feature, 
with a velocity consistent with the velocity of the companion star, is visible for 0.5 ^ 4> ^ 0.7, 
superimposed on the emission profile. 
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Fig. 11. — As in Fig. 9, for the 1996 June 10 observations, at spectroscopic phases 0.79 < <p < 0.88 
(left panel) and at phases 0.90 < < 0.95 (right panel). 
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Fig. 12. — Top panel: normalized averaged spectrum in the Ha region, obtained on 1996 June 10, 
at spectroscopic phases 0.80 < cf) < 0.92. Bottom panel: normalized averaged spectrum obtained 
on 1996 June 11, at phases 0.19 < (f) < 0.28. Wavelengths arc vacuum heliocentric. Double-peaked 
emission was observed at Ha and He I A6678. Weak [N II] A6548 and [N II] A6584 emission lines were 
probably blended with the wings of Ha. Narrow Fe I absorption lines, probably from the secondary 
star, were also seen. (They appear broadened here because the spectra have been averaged over a 
large phase interval; in particular, they appear double-peaked in the top panel because the spectra 
at ^ ~ 0.80 and (f) ~ 0.90 give the largest contribution to the average spectrum.) 



-38- 



T — I — I — r 




— 800 ' ' ' ' ^ ^ ' ' ' ' ' ' ^ ' ' ' 

0.5 1 1.5 2 

Spectroscopic Phase 

Fig. 13. — Velocities of the peaks detected in the Ha emission hne profiles, as a function of spectro- 
scopic binary phase, compared with other characteristic velocities of the system. The zeropoint in 
the velocity scale is referred to the heliocentric system; the systemic velocity of the centre of mass is 
taken as — 142.4zbl.6 km s^^ (Table 1). Open and filled circles indicate the velocity shifts of the red 
and blue peaks respectively, detected in the 1996 June spectra. Crosses indicate the peaks detected 
in the 1997 June spectra. The short-dashed line is the projected radial velocity of the primary star, 
and the dotted line is the projected radial velocity of the secondary star, as determined by Orosz 
Sz Bailyn (1997). (Smaller velocities for the two components were inferred by Phillips et al. 1999 
and Shahbaz et al. 1999). The dash-dotted lines are the projected velocities of the two edges of 
a non-circular Keplerian disk truncated at the tidal radius, as calculated by Paczyiiski (1977), for 
the system parameters determined by Orosz & Bailyn (1997); the long-dashed lines are calculated 
for the system parameters determined by Phillips et al. (1999). The velocity corresponding to the 
parameters inferred by Shahbaz et al. 1999 would be intermediate between the two. 
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Fig. 14. — Schematic elliptical disk model that can reproduce the low velocity separations of the 
Ha emission-line peaks, and the line profile variations over an orbital period as observed in 1996 
June and 1997 June. In the model, Ha would be emitted mostly from the shaded regions. We do 
see both emitting regions at all phases, but it can be noticed from the figure that at some phases 
both regions are blue-shifted (around phase 0) or red-shifted (around phase 0.5). 
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Fig. 15. — Averaged red spcctnim obtained by combining the spectra taken from 1994 August 30 to 
September 4 with the 3.9 m Anglo- Austrahan Telescope at Siding Spring Observatory. Wavelengths 
are vacuum heliocentric and the intensity is normalized to the continuum. The most prominent 
features in the spectrum are marked. 
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Fig. 16. — Averaged blue spectrum obtained by combining the spectra taken from 1994 August 
30 to September 4. Wavelengths are vacuum heliocentric and the intensity is normalized to the 
continuum. The most prominent features in the spectrum are marked. Wc notice the presence of 
broad absorption, broad (flat-topped) emission and narrow emission lines. 
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Fig. 17. — Top panel: expanded view of the Ha region in the normahzed spectrum from August 
30. Bottom panel: expanded view of the same region in the normalized, combined spectrum from 
August 31 and September 1. The narrow emission components at Ha and He I A6678 were much 
stronger on August 30 than on the following two days. 
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Fig. 18. — Ha (solid curves) and He II A4686 (daslied curves) line profiles observed between 1994 
August 30 and September 4. A red and a blue spectrum were taken in succession on each night. 
The phase indicated in each diagram refers to the observations of Ha; the phases of He II A4686 

are 4'Ha + 0.01. He II A4686 appeared generally broader than Ha, an indication that the line was 
emitted at smaller radii or closer to the disk plane. Both lines showed transitions between a single- 
peaked and a double-peaked kind of profile. The velocity zeropoint is here the systemic velocity, 
—142.4 km . On average, both lines were blueshifted with respect to the systemic velocity; see 
S6.5. 
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Fig. 19. — Top panel: medium-resolution blue spectrum obtained on 1994 September 6. Bottom 
panel: low-resolution red spectrum obtained at the same time. Wavelengths are vacuum heliocentric 
and the intensity is normalized to the continuum. The flat-topped metal emission lines and the 
broad Balmer absorption troughs observed in the previous nights had disappeared. High-ionization 
lines such as He II and Bowen N III had become much weaker; cf. Figure 16. The narrow H/3 
emission line had become stronger, but the Ha/H/? EW ratio had also increased dramatically. The 
H I Paschen lines were seen in emission. The dramatic change in the optical spectrum between 
September 4 and 6 coincided with a surge in the 20 — 100 keV X-ray flux measured by BATSE, 
and was followed by a major ejection of radio-emitting plasma on September 9. 
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Fig. 20. — Top panel: hard (20—100 keV) X-ray flux, measured by BATSE during our 1994 August- 
September observations. Bottom panel: equivalent width of the Ha, H/3 and He II A4686 emission 
lines during the same epoch. The EW of Ha on MJD 49601.36 (September 6) was 100 ± 5 A, a 
factor of 20 stronger than on the previous nights. A similar correlation between Ha emission and 
hard X-ray flux was observed in 1996 June. The EW of Ha on MJD 49590.47 has been determined 
from a spectrum taken for us on August 26 by Gary Da Costa. The values of EW on MJD 49593.47 
have been obtained from a spectrum taken for us on August 29 by Raffaella Morganti. 
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Fig. 21. — Top panel: medium-resolution blue spectrum taken with the 3.9 m Anglo-Australian 
Telescope on 1994 September 27 (HJD 2449622.867 = MJD 49622.367). Bottom panel: medium- 
resolution red spectrum from the same night (HJD 2449622.902 = MJD 49622.402). These spectra, 
obtained from the AAT archive, show that two weeks after the end of the hard X-ray flare (cf. Figure 
2 in Hjellming 1997) the system had returned to a state very similar to the one we observed before 
September 5. Cf. Figures 15 and 16 for an identification of the main features. 



